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Abstract 

Cooperative behavior is widespread in nature, even though cooperating individuals always run the 
risk to be exploited by free-riders. Population structure effectively promotes cooperation given that 
a threshold in the level of cooperation was already reached. However, the question how cooperation 
can emerge from a single mutant, which cannot rely on a benefit provided by other cooperators, is 
still puzzling. Here, we investigate this question for a well-defined but generic situation based on 
typical life-cycles of microbial populations where individuals regularly form new colonies followed by 
growth phases. We analyze two evolutionary mechanisms favoring cooperative behavior and study 
their strength depending on the inoculation size and the length of a life-cycle. In particular, we find 
that population bottlenecks followed by exponential growth phases strongly increase the survival and 
fixation probabilities of a single cooperator in a free-riding population. 


Introduction 

Cooperative behavior often provides a strong benefit for populations. But why are cooperators not under¬ 
mined by non-cooperative individuals which take the benefit but save the costs for its provision 
For higher developed organisms, there are several ways to escape this dilemma of cooperation : For individ¬ 
uals, which are capable to recognize other individuals, memorizing previous interactions and controlling 
their handling accordingly, reciprocity, and punishment can promote cooperation M- However, those 
mechanisms cannot act in organisms of modest complexity like microbes or during the early course of life 
where memory and recognition were mostly lacking. 

In such scenarios, cooperation might prevail due to the structure of the population EnDumijii; 
this idea has been studied both theoretically and experimentally in the context of kin- , group- and 
multilevel selection [IH1 El HH H21 HU HU HIJ HU HZ1 [HI SH [20l [6l [2H |22l [23l [71 12^ . If cooperators 
more likely interact with other cooperators (positive assortment), they keep most of their benefit for 
themselves and are less exploited by non-cooperators. However, due to the costs of cooperation, a fitness 
disadvantage compared to non-cooperators is still present: Positive assortment supports cooperation but 
is not necessarily sufficient to ensure its maintenance. Crucially, positive assortment can only act if 
cooperation is already established in the population such that cooperative individuals can successfully 
assort. Thus, the question remains how cooperation can emerge starting with a single cooperating mutant. 

In this manuscript, we address this issue for a generic situation of microbial populations. Cooper¬ 
ative microbes typically produce public goods whose synthesis is metabolically costly (23 HU (23 2Si 
1221 EDI EH- For example, consider the proteobacteria Pseudomonas aeruginosa and siderophore produc¬ 
tion: when iron is lacking in the environment, cooperative strains produce iron-scavenging molecules 
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(siderophores) [32] [31]. Released into the environment these molecules can efficiently bind single iron 
molecules and the resulting complex can then be taken up by surrounding bacteria. Microbial popu¬ 
lations are highly structured: several colonies form one population [331134| . New colonies arise due to 
migration into new habitats or more actively due to controlled life-cycles triggered by environmental fac¬ 
tors. For example, studies of P. aeruginosa [3511331 confirm that typical life-cycles pass through different 
steps with regularly occurring dispersal events ensued by the formation of new colonies. As the initial 
colony sizes are typically small, such dispersal events coincide with population bottlenecks. Emulating 
the dynamical colony formation, microbial cooperation has been studied experimentally by employing a 
life-cycle where new colonies are regularly formed from old ones ESI [ST]. These experiments and theo¬ 
retical work m nn El Em urn sni show that such a restructuring mechanism can cause an increase in the 
level of cooperation. 

Based on these observations, we here theoretically investigate if such a restructuring scenario cannot 
only maintain cooperation, but also allows for the evolution of cooperation from a single mutant. Assum¬ 
ing a constant population size, the onset of traits from a single mutant and its fixation have been studied 
in the frameworks of population genetics and evolutionary game theory n mi m m am m nzi mi . 
However the consequences of ecological factors such as population growth and population restructuring 
remain unclear USED]. To tackle this question for the microbial scenario introduced before, we con¬ 
sider a life-cycle model, consisting of three steps mm, see Fig. |TJ (i) In the group-formation step, 
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Figure 1. A simplified life-cycle of microbial populations. Cooperators [blue] and free-riders 
[red] in an initially well-mixed population are randomly subdivided into groups of average size uq. The 
groups then evolve separately following two main rules. First, groups with a higher fraction of 
cooperators grow faster as more public good is present in those populations. Second, in each group 
cooperators reproduce slower as free-riders as they do not have to provide the public good; for detail see 
main text. After a certain time T all groups are merged and the cycle restarts. 

a population consisting of cooperators (C) and non-cooperating free-riders (F) is randomly assorted 
into different groups (colonies), (ii) In the ensuing group-evolution step, each group evolves separately 
according to generic growth laws of microbes, (iii) In the group-merging step, colonies are merged to 
one population again. Consecutively, the cycle starts anew with the next group-formation step. The 
synchronous reformation of groups via merging all individuals into a single population follows recent 
experiments m and is obviously a simplification of natural colony formation. However, it captures the 
essence of regular occurring bottlenecks, namely rearranging colonies which have an initially small pop¬ 
ulation size. Furthermore, as it comprises a worst-case scenario for cooperation, it is suited to study the 
possible onset of cooperation in microbial populations starting with a single cooperative mutant. Starting 
with a well-mixed population with a fraction Xq of cooperators, M groups (colonies) are formed during 
the group-formation step. Successively, each group i is assigned i randomly chosen individuals of the 
well-mixed population; group sizes {r'o.i} are Poisson-distributed with mean no . The random assortment 
leads to a statistical variation in the initial fraction of cooperators, £o,i- Approximately, it is of the 
order no ■ Xo(\ — £o). We investigate this unbiased random assortment of groups as it does not assume 
individuals to be distinguishable by the sorting mechanism and again constitutes a worst-case scenario 
for cooperators. 
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After assortment into groups, each group i evolves and grows separately. The dynamics within 
groups is given by a stochastic process based on birth and death events which are characterized by the 
corresponding per capita birth rates Pg 4 and death rates ? ., where S G {C, F} denotes the trait of 
the individual in group i [SB ED- The birth rates of individuals depend on two factors, namely the 
trait of the individual and the composition of the colony the individual is living in. First, cooperators 
reproduce slower than free-riders in each colony as they have metabolic costs due to the production of 
the public good, T j ? . < T^Second, as more cooperative groups produce more public good, individuals 
in colonies containing a higher fraction of cooperators are better off, TiA < Fg • for < £j. The death 
rates incorporate the effect of limited resources and, therefore, increase with an increasing population 
size, Tg i < Tgj for i\ < Vj. For specificity, we assume the following birth and death rates which fulfill 
all conditions stated above: 


Tli =r(l +p£i)(l - 5 S , c c) and T S i = Vi/K , (1) 

where 5s,c is the Kronecker delta defined by 5c,c = 1 and 5f,c = 0. While p sets the growth advantage of 
cooperators on the colony level, the parameter c measures the metabolic costs of cooperation. The growth 
rate r = 1 is assumed to be fixed setting the time-scale of growth. The here assumed functional form of the 
growth rates reproduces the generically observed growth dynamics of microbial populations [52] : Small 
colonies grow exponentially and their size is bounded by a maximal colony size which here scales with 
K. A more detailed description of the dynamics including a discussion of the deterministic equations can 
be found in Refs. 150 ESI- Furthermore, in [ 55] the specific form of the rates (JT]) is justified by successful 
comparisons with experiments by Chuang et al. EZ3- Note also that the qualitative results presented in 
this manuscript do not depend on the specific functional forms of the growth rates, but only on the rather 
generic conditions of population bottlenecks followed by growth. 

After a regrouping time T the separated groups are merged again into one well-mixed population with 
a then changed global population size N = and a fraction of cooperators which is given by the 

weighted average, 

x = X]XI Vi ■ ( 2 ) 

i i 

The cycle then starts anew with the new fraction of cooperators, xq = x. Although the fraction of 
cooperators within each group is expected to decrease during group-evolution, an increase in the global 
fraction of cooperators is possible in principle: The disadvantage of cooperation within each group can 
be overcome by changing weights, Vi/N , in the total population. To achieve this there must be a suffi¬ 
ciently high positive correlation between group size and cooperator fraction 153 El- Such an increase of 
cooperation is an example of Simpson’s paradox mm- 

For the random assortment of groups considered here, two mechanisms promoting cooperation can be 
distinguished as previously studied |39j . First, for very small population bottlenecks, purely cooperative 
colonies might be formed where there is no conflict with free-riders ( group-fixation mechanism). Second, 
more cooperative groups grow comparably fast and thereby compensate for the selection advantage of 
free-riders (group-growth mechanism). As those mechanisms are crucial for the understanding of our 
results concerning single mutants spreading in the population, we first repeat some arguments from 
Ref. [59| and additionally introduce analytic calculations and a study of the key parameters to support 
them. The second part of this article is devoted to the main question of the paper namely whether a 
single cooperative mutant which cannot rely on benefits provided by other cooperators has the chance to 
spread in the population. 
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Results 

In the following we analyze both mechanisms in detail, starting with the group-fixation mechanism. For 
long separate evolution of groups, 1, all groups reach a stationary state: They consist of either 

cooperators or free-riders only with a group sizes of (1 + p)(l — c)K and K , respectively. The global 
fraction of cooperators is then 


(1 +p)(l - c)P c 

(l+p)(l-c)P c + (l-P C ) ' 


(3) 


Pc denotes the probability for a group to consist of only cooperators after assortment. In first order, only 
initially purely cooperative groups contribute to Pc while all mixed groups are taken over by free-riders, 
such that 


1 _ _ Vi 1 

p c = — 1 Y.^ + °( 1 <) = 

u,=l 1 


e n 0 x 0 _ 1 

e n ° - 1 


+ 0 ( 1 ). 


If x' exceeds the initial fraction of cooperators, xo, the group-fixation mechanism is strong enough to 
overcome the advantage of free-riders. As Pc increases with the initial fraction of cooperators, there is 
an unstable fixed point x* u , implicitly defined by x* u = x' = Xq in Eq. For initial fractions, Xq , above 
x* , a purely cooperative population is reached after several regrouping events. In contrast, when starting 
below, Xq < x*, cooperators become extinct in the population. This bistable behavior is illustrated in 
Fig- i where depending on the initial value Xq the global fraction of cooperators is shown after a large 
regrouping time T = 20. 

In contrast to the group-fixation mechanism, the group-growth mechanism acts for small times, where 
groups strongly grow. As cooperation enhances the growth speed of colonies, more cooperative groups 
have a larger weight in the average ([2]) even though < 0 holds in each group. Depending on the param¬ 
eters, this positive effect is able to compensate for the selection disadvantage of cooperators. This can be 
quantified performing a van Kampen expansion of the master equation, see supplementary information. 
For binomial distributed groups, the change in the fraction of cooperators at time t = 0 is given by, 


d 

dt 


x oc 


[—c(l +px) + p/n 0 \ x(l 


x). 


(4) 


The first term accounts for the selection advantage of free-riders for the growth advantage of more 
cooperative groups while the second one reflects the growth advantage of more cooperative groups. The 
initial change, -^x, is larger if xq is small meaning that higher selection disadvantages can be overcome. 
Under regrouping after time T, the group-growth mechanisms results in a stable fixed point x* s . As 
the group-growth mechanisms relies only on variance in group composition, but not on the existence of 
purely cooperative groups, it acts for much stronger population bottlenecks, no, than the group-fixation 
mechanism does. However, as it is caused by population growth it can only act for short regrouping times 
T. 

Repeated regrouping corresponds to an iterative map. The underlying dynamics (group formation 
and group evolution) result in an effective drift: Ax = x{T) — Xq. Depending on the strength of both 
mechanisms, and thus no and T, five distinct fixed point scenarios can be distinguished; examples for 
the corresponding stability plots (Ax) are shown in Fig. [2j Besides the discussed bistable [only group- 
fixation] and stable coexistence scenario [only group-growth], there can be a bistable coexistence scenario 
[group-growth and fixation mechanism]. In addition, there are the scenarios of only cooperation and only 
defection , where the sole stable fixed points are x* = 1 or x* = 0, respectively. 

To this point we have seen that the internal dynamics on the intra-group level result in an effective 
drift for the iterative map, which can support cooperation. We now want to examine whether this drift 
allows for the evolution of cooperation. In other words, we investigate if a single cooperative mutant 
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Figure 2. Parameter dependence of the different cooperation scenarios. Depending on no 
and T five different scenarios with different fixed point behaviors arise. For each scenario, we show an 
exemplary drift diagram on the right, where the change in the fraction of cooperators after one life-cycle 
is shown depending on the initial fraction of cooperators, Ax(xq)- Dots correspond to simulation results 
of the transition points and the lines are guides to the eye to separate the different scenarios: Pure 
cooperation with a stable fixed point at 1, i.e. for a purely cooperative population (no = 4, T = 1.5), 
coexistence with a stable fixed point at 0 < x* < 1 (n o = 6, T = 1.8), bistability with an unstable fixed 
point at 0 < x* < 1 (no = 5, T = 20), bistable coexistence where both an stable and an unstable fixed 
point are present (no = 4, T = 5.5) and only free-riders. The black square corresponds to the 
parameters studied in Fig. [3] ( no = 5, T = 4). Other parameters are p = 10, K = 100, and c = 0.1. 


can survive and spread in a free-riding population. In Fig. [3] we show two exemplary trajectories for 
the time evolution of the fraction of cooperators. The parameters T = 4 and n o = 5 corresponds to the 
coexistence regime in Fig.[2](black square). The simulations are performed by placing a single cooperative 
mutant in one of the M = 200 groups; this corresponds to an average initial fraction of cooperators of 
xq = l/(noM) = 0.001. In Fig. [3] A, a coexistence trajectory is shown while in Fig. [3] B cooperators go 
extinct at the fifth regrouping event (t = 20). There are strong differences between these trajectories and 
the ones known from non-linear dynamics models without an iterative map which are commonly used to 
analyze coexistence scenarios: Due to the dynamics during group-evolution, the fraction of cooperators 
does not only fluctuate slightly at the fixed point, but oscillates around it. The reason is the interplay of 
the group-growth mechanism with the ensuing decline in the fraction of cooperators due to the selection 
disadvantage of cooperators in each group: After each regrouping event, the group growth mechanism 
causes an increase in the level of cooperators. Later on, when population growth slows down as groups 
are already close to their carrying capacity, the group-growth mechanism is not effective anymore. Thus, 
the fraction of cooperators declines due to their growth disadvantage within each group. This dynamic 
pattern is then repeated in each life-cycle causing the oscillatory behavior. But not only oscillations 
increase the variability, also demographic fluctuations are particularly strong. The reason is that small 
variations in the fraction of cooperators are exponentially amplified during group growth. Therefore, 
for large times (t > 100 in Fig. [3] A) when the steady state is already reached the exact position of 
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Figure 3. Two exemplary realizations for the evolution of cooperation from a single 
mutant. The parameters, T = 4 and no = 5, correspond to a point in the coexistence regime indicated 
by the black square in Fig. [2j In panel A a coexistence trajectory is shown, while panel B shows 
extinction of cooperators. For small times (zoom in inset panel A, and panel B) cooperators have to 
survive a highly stochastic process: Cooperators start with only one mutant, the population is subject 
to exponential growth which amplifies fluctuations, and regrouping events add an additional level of 
stochasticity. In addition, also for larger times the trajectory oscillates around the fixed point (for a 
detailed explanation see main text) and is therefore prone to extinction. 


the maxima and minima in each life-cycle substantially varies between different regrouping events. In 
addition, also the regrouping mechanism increases the randomness, i.e due to the stochastic reformation 
of groups the fraction of cooperators can change during this process. In the inset of Fig. 0A the first eight 
regrouping steps are magnified. For instance during the second (t = 8) and fourth (t = 16) regrouping 
events a substantial drop in the fraction of cooperators is present, while for other regrouping events it is 
hardly changed. 

This analysis already shows, that the question whether a single mutant can survive is much less 
trivial than suggested by the phase diagram [Fig. [2]: Mutants have to overcome two hampering factors 
which diminish the positive drift of the iterative map. First, in each group free-riders are favored over 
cooperators due to the costs to produce the public goods c: During group-evolution, the fraction of 
cooperators declines in each initially mixed group. Thus if group-evolution was not interrupted by the 
formation of new groups, cooperators would die out in those mixed groups. Therefore, in all parameter 
regimes (including the ones with a stable cooperation fixed point) the fixation probability of a free¬ 
rider is higher than the one of a cooperator during group-evolution. Second, even if cooperators have 
survived group-evolution, it is not assured that they survive regrouping as this process follows a random 
distribution. Especially for small fractions of cooperators, cooperators die out more likely than free-riders 
during regrouping. Taken together, it is not obvious whether the advantages of cooperators are sufficient 
to enable the evolution of cooperation from a single mutant. 

In the following we study this issue by analyzing the fixation probability Pf, and the survival probability 
P s . Pf = -Pfix.c denotes the probability that a single cooperative mutant will take over the entire popu¬ 
lation, whereas P s = 1 — Pr x ,f is the probability that free-riders have not fixated yet in the population. 
Both probabilities are time-dependent, but reach quasi-stationary values under repeated regrouping, see 
supplementary information. Again starting with one cooperator, we performed stochastic simulations of 
the population dynamics, each realization with 200 regrouping events. In Fig. [4j the probabilities P s 
(panel A) and Pf (panel B) depending on ng are shown for different regrouping times T. Due to the 
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Figure 4. Survival (panel A) and fixation probability (panel B) of cooperation. When 
starting with a single mutant, the survival probability of cooperators, P s , and their fixation probability, 
Pf, are shown after 500 regrouping events. Interestingly both probabilities are fairly large over a large 
parameter regime. Different colors correspond to different regrouping times: 

T = {2 (*),4 (*),5 (*),6 (•), 10 (•)}, M = 200, number of realizations R > 10 5 . Other parameters: 
p = 10, K = 100, c = 0.1. Black lines are calculated with our analytic approximation, Eq. <© for T = 2 
(solid) and T -» oo (dashed). 


life-cycle dynamics the fixation behavior here is more intricate than for standard evolutionary dynamics. 
In those standard models, it is the population size only that determines the weight of fluctuations with 
respect to the deterministic drift [2 . In the limit of large population sizes fluctuations become irrelevant 
and a stable fixed point is reached with absolute certainty. In contrast, for our life-cycle model, the stabil¬ 
ity diagram and the deterministic drift themselves depend on the initial population size n$: cooperators 
are less favored for higher values of no. For small no, the regime of pure cooperation with a stable fixed 
point at x* = 1 is present, cf. Fig. [2j In accordance with that both fixation and survival probability 
are equal and fairly large. With increasing initial population sizes Uq the strength of the group-growth 
mechanism is reduced and we observe a decline in the survival and fixation probabilities. This decline 
is faster for larger regrouping times since cooperators stay longer in the evolution step and the thereby 
created selection disadvantage, s, accumulates during the time interval [0, T]. Further increasing initial 
groups sizes cause a steep drop in the fixation probability. The reason is that upon increasing no one is 
leaving the regime of pure cooperation, cf. Fig. [2] red line. Hence, the fixation probability of cooperators 
dramatically declines. While the so far discussed fixation probabilities only depend quantitatively on 
the regrouping time T, the survival probabilities are also qualitatively influenced by this parameter. For 
small T, the system is in the coexistence regime which implies that cooperators neither fully take over 
the population nor die out easily. Thus the survival probability shows only a slow decay with no in this 
regime. In contrast, for larger regrouping times, the dynamics becomes unstable. Because the dynamics 
start with a single cooperative mutant, the group-fixation mechanism is not efficient as it is unlikely that 
the threshold above which cooperators are favored is reached. Therefore, the survival probability declines 
rapidly with increasing no- 

Let us finally discuss how the survival probabilities can be understood, based on analytic arguments. 
To this end, we approximate the dynamics by considering the first life-cycle only, i.e. the first group- 
evolution steps and the ensuing formation of new colonies. As most extinction events of cooperators 
happen at the beginning, this approximation captures the extinction dynamics qualitatively correctly, 
see supplementary information. The success of a cooperator crucially depends on the size of the group 
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Figure 5. Robustness of the survival probabilities for different model parameters. Panel A 
shows the impact of exponentially distributed regrouping times depending on tlq. In panel B the impact 
of varying selection disadvantages is investigated for p = 10 and no = 5. In panel C the survival 
probability for different values of p scaling the growth advantage mediated by the public good is shown. 
Even for small values of p the group-growth mechanism is effective and support the onset of 
cooperation. If not varied, no = 5 and c = 0.1. Other parameters are as in Fig. 4. 


it is living in. For larger groups it has to compete with more free-riders and its survival chances are 
diminished. As the group sizes are Poisson-distributed, the probability that a cooperator emerges in a 
group of size k is given by n,Q~ 1 /(k — l)\e ~ n °. The probability that cooperators survive the first regrouping 
step is the probability that not all Muq newly formed groups are purely free-riding. For a realization with 
a fraction of x cooperators before regrouping, this probability is given by 1 — e~ xn ° M . Taken together 
this leads to the overall survival probability 


P s = e - no J 2 n o~ 1 /(k-iy-(l-e~ x{k:r)noM ) (5) 

fc=o 

where (..) denotes an average over all possible realizations. For large regrouping times, cooperators die 
out in all initially mixed groups, i.e. x = 0 for k > 1. Thus, the survival probability simplifies to e ~ n °, 
which is in perfect agreement with our simulation results; see Fig |4] (dashed black line). For smaller T, 
the group-growth mechanism is dominant, and cooperators also survive in initially mixed groups leading 
to even higher values of the survival probabilities. For T = 2, where the group-growth mechanism results 
in a stable fixed point, we compare the approximation with the simulated survival probabilities and find 
good agreement; cf. Fig [4] (solid black line). 

To further support the idea that ecological factors can explain the onset of cooperation, we study 
the robustness of our findings against changes of the model parameters. In particular we vary p scaling 
the benefit of the public good, and c measuring the metabolic costs due to the production of a public 
good. In Fig. [5j we show the survival probabilities for different regrouping times and n o = 5 depending 
on c (panel B) and p (panel C). Interestingly, already for comparably small values of p the positive 
impact of a public good on the population size notably supports cooperation and thereby increases the 
survival probability. In contrast increasing costs c decrease the survival chances, but the benefits of 
cooperation are strong enough to compensate for selection disadvantages of up to 20%. In addition, we 
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also demonstrate that our results do not depend on the assumption that the regrouping time is fixed. In 
Fig. | A, survival probabilities are shown for regrouping times exponentially distributed with T according 
to p{T) = i exp[— T/T\. Again, the survival probabilities are fairly large and our conclusion that a simple 
life-cycle supports cooperators strong enough to overcome their selection disadvantages is not changed. 
All in all, we therefore conclude that our results are robust against changes in the details of the parameters 
and model assumptions and thereby constitute a possible mechanisms enabling the onset of cooperation. 


Conclusion 

In this article, we studied the impact of ecological factors like population growth and population bottle¬ 
necks on the evolutionary dynamics of cooperating individuals. Our main findings can be subdivided into 
two blocks: First we analyzed the evolutionary dynamics acting on cooperators already abundant in a 
population, second we studied the survival chances of single mutants emerging in a purely non-cooperating 
environment. 

In the first part, we were employing our recently introduced model [59] to study how a restructuring 
mechanism combined with typical growth conditions influences the evolutionary dynamics of public good 
producing bacteria. Depending on the inoculation size no and the regrouping time T, regimes of stable 
cooperation, coexistence between cooperators and free-riders and bistability emerge. Those regimes 
arise over a broad parameter regime even though the worst case scenarios for cooperators are assumed 
whenever model assumptions have to be made. Therefore, we believe that the mechanisms still apply 
in more realistic evolutionary scenarios where for example reassortment is not completely random or 
public goods are not equally distributed between all individuals [9j. Other studies focussing on different 
aspects of the interplay between evolution and ecological factors support our findings: For instance, 
the impact of exponential growth following bottlenecks in infinite populations was shown to support 
cooperators [19] . In Ref. [50] the competition between groups sharing a bounded global populations 
size and thereby competing for resources was investigated, while in Ref. m the impact of mutations 
on the beneficial effects of population growth for groups starting with only one individual was studied. 
In addition, also the frequency of ecological disturbance and resource supply plays a crucial role for the 
resulting level of cooperation [53, [55]. Furthermore, beneficial effects for cooperators were also found 
when reassortment is not random but environmental driven I5EI E3 Eg. All studies emphasize the 
important role of ecological factors for understanding cooperation. In contrast to the mentioned studies, 
we focussed on a description of bacterial growth, starting with an exponential growth phase reaching 
a carrying capacity later on. The different growth regimes (exponential growth and stationary state) 
influence the evolutionary dynamics differently: Both related mechanisms (group-growth and group- 
fixation) favor cooperators, but as confirmed by analytic arguments presented above, the strength of 
the group-growth mechanisms increases with small .To while the strength of group-fixation mechanism 
decreases. Therefore, both mechanisms can be associated with two different fixed point scenarios (stable 
and unstable). We present a full parameter study of the ensuing regimes of cooperation (fully cooperative, 
coexistence, bistability, bistable coexistence and purely free-riding) for both key parameters, rig and T. 

In the second part of our paper we focussed on the question whether those beneficial mechanisms can 
also explain the onset of cooperation from a single mutant. Similar questions were extensively studied 
for evolutionary dynamics without regrouping where the factors influencing fixation probabilities for 
neutral, beneficial and deleterious mutations were investigated carefully, see e.g. 59] 160] |H]. However, 
the non-iterative map caused by regrouping makes a new study essential as it alters many aspects of the 
evolutionary dynamics. Here, we demonstrate that ecological factors increase the survival and fixation 
probabilities of cooperators substantially. In particular, the group-growth mechanism allows for the robust 
establishment of cooperation as it does not rely on a threshold fraction of cooperators to act effectively. 
Remarkably, the probability for a single mutant to succeed decreases only slowly with increasing n o, and 
growth thus allows the onset of cooperation without the requirement of extremely narrow population 
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bottlenecks. Due to the robustness against parameter changes and the worst case assumptions employed, 
our model provides a proof of principle that ecological factors might explain the onset of cooperative 
behavior. Once cooperation is established in a population, more advanced mechanisms, which rely on 
cooperators already present in a population, like kin discrimination or other active forms of positive 
assortment, may evolve to further stabilize cooperative behavior, see e.g. [II EH E31 El ESS ESI HU- 
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Supporting Information 

In this supplementary material, we show how van Kampen’s system size expansion can be used to derive 
the initial slope of the fraction of cooperators, Eq. (4) main text. In addition, we present data which 
confirms the validity of Eq. (4) and approximations made for the calculations shown in the main text. 


1 Van Kampen Expansion 


As already described in the main text, the group growth mechanism is a fluctuation driven effect. To 
describe it, we perform an Omega expansion [T] of the master equation (ME) defined by the birth and 
death rates Eqs. (1) (main text). The results obtained in this section hold for a fixed initial group 
size, where the random assortment only influences the composition but not the number of individual 
of a colony. In most of our simulations (all except Fig. SI), we use Poisson distributed group sizes 
and fluctuations in the initial composition of a group are additionally enhanced by variation in group- 
size, see main text. This full solution can be calculated by weighting the result for fixed colony sizes 
by a Poisson distribution with mean, no- But as the important parameter dependences are already 
present when considering only binomial distributed groups and the expression are less lengthy, we focus 
on this simpler case in the following. Before we start with the calculations, we also want to make 
another comment. Normally, the condition for successfully performing the Omega expansion, is a stable 
deterministic attractor. However, here we do not want to analyze the resulting Fokker-Planck equation 
in detail, but only use it to determine an estimate of the strength of the group-growth mechanism at 
t —>• 0 (if the mechanism is not strong enough to overcome the selection pressure at the very beginning, 
it will also not happen later as fluctuations become smaller with increasing time). Therefore, the strict 
assumption of a stable fixed point can be relaxed, which leads to good results as we demonstrate at the 
end of this section. The ME describing the evolutionary dynamics is given by, 

d t P(N Cl N F ) =[(Ec - l)r(l +P Nc N ° Nf ) n c + (E^ - lHl+p ^^ Kl + s)N F 

+(E+-1 )N c Nc + Nf + (E+ -1)N F N ° + Nf }P(N g N f ), (SI) 

where P(Nc, N F ) is the probability to find a group containing Nc cooperators and N F free-riders. The 
ladder operators Eg are defined by Eg/(TVs) = f(Ns ± 1). 

Now, within the van Kampen system size expansion, the stochastic variables, Nq and N F , are sub¬ 
stituted by a linear combination of functions describing their rescaled means, c(t) and f(t), and new 
stochastic variables, £ and /x, accounting for fluctuations, 

N c =Kc(t) + VK£, 

N f =Kf(t ) + y/K/i. (S2) 


In this expressions fluctuations are already weighted according the well-known square root dependence 
on the population size [5] which scales with the carrying capacity K. 

Employing Eqs. ( |S2| ) in Eq. (SI) and expanding the resulting expression in the carrying capacity I\ 
leads to the mean-field equations for c(t) and f(t) in order \J~K and a Fokker-Planck equation (FPE) for 
P(£, n) in higher orders. As the FPE is lengthy and does not contributed to the understanding of our 
calculation, we do no write it down explicitly here. In Table [l] we summarize all changes arising with 
the substitution and the expansions which are used to transform the ME to a FPE equation. 

Next, we analyze the mean values and second moments of the fluctuations. They can be calculated 
by multiplying the FPE with the respective fluctuation variable and integrating over both, f d^dfj, £... 
with £ = {£, u}. At time t = 0, the means of the fluctuations vanish according to the initial conditions 
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original form 

substitution and expansion 

d t P(N c ,N F ) 

dtP{£, ») - VK(c(t)dsP(Z, p) + /(i)^P(£, p) 

Ef 

11,2 ~7K d i + 

N s 

Ks{t) + VK( 

N s 

K 2 a{t ) 2 + 2 K 3 / 2 a{t)C + K ( 2 

N c 

Nq-\-Nf 

X W + m(t)VK K 1 * 


Table 1. Table summarizing the substitutions and expansions which are employed to transform the 
ME into a Fokker-Planck equation. The stochastic variable ( represents the fluctuations £ = {£,/r}. 
The means are either given by a(t) = {c(f), /(£)} or x(t) = and m{t) = c(t) + f(t). 


log (t c ) 



Figure SI. Condition for the increase in the level of cooperation caused by the group-growth 
mechanism. The cooperation time is shown for different values for s and no- The black line corresponds 
to the condition, (S7). The other parameters are p = 10, r = 1 and K = 100. 
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Figure S2. Temporal change of the survival probability of cooperators, Ps, when starting with a single 
mutant. Change occurs only during the first regrouping events. Lines show Ps for different initial group 
sizes no and fixed regrouping time T = 2. 
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(Co) = {p o) = 0. Thus the time evolution of the means is determined by only the second moments and 
given by, 


(C) — H- 7 == [(1 — a ; o) 2 (Co)“2a;o(l — ®o)(C/ i o)+a ; o(/ x o)] > 

movA 

(A) = _ rp ^ ~tf) [(l-xo) 2 (Co)-2xo(l-a;o)(CMo)+®o(Mo)] , 

moviv 

(S3) 


where we used mp = 
binomial distribution, 
moments, (£g) = (mo) 


c(0) + /(0) and Xq = Cp/rrip. As the composition of each group is drawn from a 
the variance of this distribution sets the initial condition for the second fluctuation 
= — (Cmo) = fix o(l — Xq). Employing this in Eq. (S3) leads to, 


, a _ rpx 0 (l - x 0 ) 

w Vk ’ 

(S4) 

v A 

The reason for the different signs of both differential equations is that the global growth rate 1 + px 
is positively correlated with the fraction of cooperators x while it is negatively correlated with the 
fraction of free-riders, 1 — x. Hence, cooperators are favored by fluctuations. To compare the strength of 
this fluctuation driven mechanism with the deterministic selection disadvantage, we analyze the average 
fraction of cooperators. To first order in 1/A, it is given by, 


(x) =i 


x+ 


71+ -±=(£ + n) 


* x (t) + ((1 - x )(0 ~ x(n)) ■ 


(S5) 


Differentiating this expression, neglecting higher order of s and 1/A, and using Eqs. (S4) leads to, 


^l( x ) = - sr(l +px 0 )x 0 (l - x 0 ) + rpXo ^ x °"> = 

at 777-0 A 

= - sr(l +px 0 )x 0 (l - x 0 ) + rpXo<yl —^1. (S6) 

n 0 

To confirm the validity of this calculation, we determine the parameter regime, in which the positive 
second term can overcome the selection disadvantage caused by cooperators ^(x) > 0. This condition 
holds if, 


s > 


n 0 (l +pxo) 


(S7) 


In Fig. |S1[ we compare this condition to simulation data, for varying s and no- We use the cooperation 
time, the time until the fraction of cooperators drops under its initial value, as a measure. If the 
cooperation time is zero, the advantage due to fluctuations is not large enough to overcome the selection 
disadvantage, while for finite values of the cooperation time, the group-growth mechanism is sufficient 
to compensate for slower growth of cooperators. As shown in Fig. SI Eq. (S7) is in agreement with 


the position where the cooperation time vanishes. Note that the analytical solution derived before gives 
a lower bound for the the fluctuation driven enhancement of cooperation: We only take into account 
the initial fluctuations given by a binomial distribution. In addition, demographic fluctuations further 
enhance cooperation, explaining the systematic shift in Fig. SI. 
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2 Time-Dependence Survival Probability 

In Fig. 3 in the main text, we show the extinction probabilities after 100 regrouping steps. Here, we want 
to confirm that after this time, already a quasi stationary state is reached. In Fig. |S2| the time evolution 
of the survival probability of cooperators is shown for exemplary parameters. The data clearly shows 
that, as already mentioned in the body of this paper, the first regrouping steps mainly determine the 
fate of the population. If a single cooperative mutant survives those first grouping events, its extinction 
probability almost vanishes. Besides validating that after 100 regrouping steps already a quasi steady 
state is reached, the presented data also justifies our analytical approach to estimate survival probabilities 
based only on the first regrouping event. 
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